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1. Introduction 

The classification of tlie universality classes of non-equilibrium phase transitions is still 
an open problem of statistical physics [H [21 El H] ■ In equilibrium conformal invariance 
(CI) [HI m [7] enables this in two dimensional critical systems as the consequence of 
a larger group (the CI group) than the mere scale transformations. Recently the 
generalization of the generators of CI (albeit without invariance under time-translations) 
are proposed for anisotropic, dynamical models [3, [9]. The corresponding invariance is 
the so-called local scale-invariance (LSI). Since it is supposed to be the extension of the 
dynamical scale transformations for such systems it may serve as a convenient tool for 
classifying universality classes of non-equilibrium systems as well. 

The quantities of main interest are the two-time autocorrelation function C{t,s) 
and the auto-response function R[t,s), which describe ageing phenomena (for recent 
reviews see [10]) 

C{t,s)= (0(t,f)0(s,f)) (1) 
5(0(t,r)) 



R{t,s) 



6h{s, r) 



= (0(t,f)0(s,f)) (2) 

h=0 

where and are the fields in the Janssen-de Dominicis formalism [Til [12] and h is the 

magnetic field conjugate to 0. For t,s ^ oo and y = t/s > 1 one expects the scaling 
forms 

C{t,s) = s-''fc{t/s) (3) 

R{t,s) = s~'-yRit/s), (4) 

where a and b are ageing exponents and fc and /r are scahng functions such that 
fc,Riy) ^ y~^c,R/z for y 3> 1. Here Ac and \r are the auto-correlation [13] and auto- 
response [H] exponents respectively and independent of equilibrium exponents and the 
dynamical exponent Z (defined as usual Z = i'\\/i'±). 

As in case of CI one expects that LSI fully determines the functional form of the 
scaling functions. Henkel et al. derived R{t,s) in general and the form of C{t,s) for 
Z = 2 by identifying the quasi-primary operators of the theory [T^ITB] . The generalized 
form of R{t, s) takes into account the difference between physical observable defined in 
lattice models and the associated quasi-primary scaling operators of the underlying field 
theory as well. This ansatz looks as 

/^x I+o'-Xr/Z X -l-a' 

i?(t,.) = .-^-m , (5) 



.8/ \S 

where a' 7^ a is an independent ageing exponent in general. Some systems with 
detailed balance symmetry has been analyzed recently and found to satisfy ([5|) 
[9l[T8l[T9ll20l[2lll22l[23] with a 7^ a' . On the other hand renormalization-group results for 
some important universality classes concluded that a = a' should be hold. In particular 
explicit two-loop field-theoretical computation of i?(t, s) for the 0{N) universality class 
and Model A dynamics at the critical point claim a = a' 
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Recently numerical simulations of the non-equilibrium contact process (CP) did 
not satisfy that form completely [26] and Hinrichsen argued that LSI is not a generic 
property of ageing phenomena but is restricted to diffusive models {Z = 2) or above 
the upper critical dimension. In [16] Henkel et al. suggested that there is crossover in 
case of nonequilibrium critical dynamics because both the ageing regime (t — s ~ 0(s)) 
and the quasi-stationary regime {t — s << s) display scaling with the same length scale 
L{t) oc t^/^. For a more detailed discussion of these results see a very recent review |17j . 

In this paper I present simulation results for an other nonequilibrium critical model, 
the parity conserving (PC) nonequilibrium Ising model (NEKIM) in 1 + 1 dimensions. I 
provide numerical evidence that in this model C{t,s) and R{t,s) can be fitted with 
the forms Eqs.QjQ hence this nonequilibrium critical model exhibits LSI scaling 
invariance. 

2. The PC class NEKIM model 

The NEKIM model has been introduced and analyzed first by Menyhard [27] as a 
generalization of the Kinetic Ising model [28] by adding spin-exchange updates in 
between the spin-flipp steps of the Glauber Ising model. In one dimension the domain 
walls (kinks) between up and down regions can be considered as particles. The spin-flip 
dynamics can be mapped onto particle movement 

Tii-TTi ~ .o^o. (6) 

or to the annihilation of neighboring particles 

Tit^TTT ~ oo->«« (7) 

Therefore the T = Glauber dynamics is equivalent to the annihilating random walk 
(ARW). This is a double degenerate phase, an initial state decays algebraically to the 
stationary state, which is one of the absorbing ones (all spins up or all spins down, 
provided the initial state has an even number of kinks). By mapping the spin-exchange 
dynamics in the same way more complicated particle dynamics emerges, for example: 

TTii'^nn ~ o.o^... (8) 

one particle may give birth of two others or three particles may coagulate to one. 
Therefore this model is equivalent to branching and annihilating random walks with 
even number of offsprings [301 El] ■ By increasing the spin-exchange a second order phase 
transition takes place [27] for the kinks from absorbing to active state, which belongs to 
the parity conserving (PC) universality class [29 l [3T | [32]. 

In [36] this model has been investigated by high precision cluster simulations with 



the parameterization 

Pe. = i- 2r (9) 

Wi = T{l-6)/2 (10) 

Wo = T{l + 5) (11) 
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originating from the Glauber Ising model [2S] . In this paper I present simulations at the 
critical point determined in previous works [271 ESI EH ESI EE] . The parameters chosen 
are: F = 0.35, pex = 0.3, 6c = —0.3928(2). Here the kink (n^ G (0, 1)) density decays as 

< n,{t) >oc t"°-285(2) (12) 

as can be seen on the inset of Fig{Tl In previous works the NEKIM algorithm introduced 
as follows. The spin-flip part was applied using two-sub-lattice updating. Following this 
states of the spins are stored and L (L is the size of the system) random attempts of 
spin-exchanges are done using the stored situation of states of the spins before updating 
the whole lattice. All these together was counted as one Monte Carlo time-step (MCS) 
of updating (throughout the paper time is measured by MCS). 



3. Simulations 

Time-dependent simulations were performed in L = 2 x 10^ — 10^ sized systems with 
periodic boundary conditions. The runs were started from fully ordered kink state 
(nj(0) = 1) i.e. alternating up-down spin configuration (sj G (—1,1)). I followed the 
quench towards the critical state and measured the kink (order parameter) density 

< ni{t) >oc r" , (13) 

the kink-kink autocorrelation 

C{t,s) =<niit)n,{s) > , (14) 

and the auto-response function, by flipping a spin at random site / at time s generating 
a kink pair out of the vacuum 

R{t, s) =< ni{t)ni+i > - < n-(t)?2^+i > \,>^s):=-s,{s) (15) 

The simulations were run for several values of waiting times s = 256, 512, 1024, 2048, 4096 
and the scaled autocorrelation C{t,s)t~'^°' is plotted on Fig. [T]with the assumption of 
the form Eq.([3]). Good data collapse (within error margin of the simulations) in case of 
systems of sizes L = 2x 10"^ could be achieved for the whole region, however for larger s 
and t values small deviations form the collapse could also be observed. By investigating 
larger systems this proved to be finite size effect. The curve on Fig. [1] for s = 4096 
shows the result of L = 10^ simulations. In the asymptotic t/s — * oo limit it can be 
fitted by power-law, corresponding to the density decay of the PC class [36]. This 

suggests the scaling exponents: b = 0.570(4), Xc/Z = 0.285(2). By inserting the value 
of the dynamical exponent of the PC class Z = 1.75(1) [3j on obtains Xc = 0.498(2). 
This exponent agrees with that of the autocorrelation exponent of spins A = 1.50(2) 
[35] for (t/s) oo 

A{t,s) = < S^{s)s^{t) > 

= f{t/s) oc (t/.)-(^-'^+i-V2)/z (16) 

since rj = 1.01(1) and a hyperscaling-law connecting time dependent spin and kink 
exponents at the PC transition point derived in [3l]. Fitting for the connected 
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Figure 1. Autocorrelation C{t,s) of the critical NEKIM model for several values of 
waiting time s as a function of the scaling variable t/s. The dashed line has the slope 

-0.285 a. The inset shows the kink density decay p{tf-^^^ up to t = 409600 MCS 

in a system of size L = 10^. 



autocorrelation function defined as r(t, s) = C(t,s) — N{t)N{s) resulted in Xq/Z = 
1.9(1) 

The auto-response function has been found to exhibit similar nice data collapse 
by plotting R{t,s)t^'^'^ as the function of y = t/s (Figj2]). However in [26] Hinrichsen 
discovered that in case of the CP model deviations from the LSI scaling form of R{t, s) 
([5]) may occur for y ^ 1. To see this I plotted R{t, s)s^'^'^y^{y — 1)^ suggested in [26] 
as the function of y — 1. Now one may see agreement with Eq.(l5]) if the curves fitted 
with the parameters collapse and are horizontal for all y values. The best agreement 
could be achieved with A = —1.3(1) and B = —0.57(1) plotted in the inset of Figl2j 
By increasing s the observed deviations from LSI scaling for y —>■ 1 occur at smaller y 
values, suggesting corrections due to the microscopic reference time s. This is different 
from the case of the CP, where all such curves collapsed. Assuming the general form 
Eq-dS]) fitting results in: a = a' = —0.430(2), Xr/Z = 1.9(1) dynamical exponents with 
a validity of more than three decades. 
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Figure 2. Auto-response R{t, s) of the critical NEKIM model for several values of 
waiting time s as a function of the scaling variable t/s. The dashed line is a fit with 
the general form ([5|). In the inset the rescaled auto-response function is plotted in 
such a way that possible deviations from the LSI scaling are more easily observable for 



4. Conclusions 

In conclusion numerical simulations of the parity conserving NEKIM model in Id 
support local scale invariance at the critical point. In contrast to the contact process 
(belonging to the directed percolation class [Ij) corrections to scaling due to the 
microscopic reference time vanish in the s ^ oo limit. Both the autocorrelation and 
the auto-response functions can well be described by the functional forms of LSI, only 
negligible dependence on the system sizes has been detected within error margin of the 
numerical simulations. Further sources of deviations may come from the value of a and 
the location of the critical point. The same analysis done for —6c = 0.3925, 0.392, 0.391 
and using a = 0.286, 0.287 did not result in viewable deviations in the figures and 
the fitting parameters. The auto-response function scales in such a way that a = a'. 
Numerical estimates for the Xc,g,r exponents are determined and Xq = Xr = Z{l+a)+d 
in agreement with the scahng hypothesis of [37]. The results support the conjecture of 
Henkel [15] that LSI can be extended to other nonequilibrium critical systems besides 
diffusive models {Z = 2) and models above the upper critical dimension. 
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